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Abstract 

We present model independent constraints on the masses and couplings to 
fermions of B and L conserving leptoquarks. Such vector or scalar particles could 
have masses below 100 GeV and be produced at HERA; we list the generation 
dependent bounds that can be calculated from rare lepton and meson decays, 
meson-anti-meson mixing and various electroweak tests. 
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1 Introduction 



A leptoquark is a scalar or vector particle that couples to a lepton and a quark. It 
may or may not have well-defined baryon and lepton number, depending on the choices 
of coupling. Those with B violating interactions would in general mediate proton decay, 
so their masses are expected to be very large (~ 10 15 GeV). In this paper we only 
consider leptoquarks with B and L conserving renormalizable couplings consistent with 
the symmetries of the Standard Model. 

There are no interactions involving a quark, a lepton and a boson in the Standard 
Model. There is a scalar Higgs doublet with electroweak quantum numbers, and vector 
bosons that are either coloured or charged, but no boson carrying colour and charge. 
This is a reflection of the fact that classically the leptons and quarks appear to be 
independant unrelated ingredients in the Standard Model. However, in each generation 
of the quantum theory, they have equal and opposite contributions to the hypercharge 
anomaly, which must vanish for the quantum theory to make sense. It would therefore 
seem natural to have interactions between the quarks and leptons in any extension of 
the Standard Model, and, in consequence, bosons coupling to a lepton and a quark. 

At ordinary particle-anti-particle colliders, leptoquarks could be pair-produced, lead- 
ing to a lower bound on their mass of roughly half the accelerator centre of mass energy. 
However, heavier ones can only be exchanged in the t channel, giving a small cross- 
section for the interactions they mediate. This problem is not present at a machine 
colliding electrons and protons, since a quark and an electron could form an s channel 
leptoquark; HERA, the e-p collider at DESY, which has already published some bounds 
on leptoquarks f48|, should be able to see |], || those with rrii q < 314 GeV (the col- 



lider centre-of-mass energy) as peaks in the x distributions of inclusive processes, if such 
leptoquarks exist. They also expect to see virtual effects due to heavier leptoquarks 

Present and upcoming results that could give constraints on leptoquarks were dis- 
cussed in J5J. At low energies, leptoquarks could induce two-lepton two-quark interac- 
tions, like those mediated by the electroweak four-fermion vertices. This suggests that 
the leptoquark yukawa coupling squared (= A 2 ), divided by the mass squared (= mf q ) 
is at least as small as Gp. However, for certain combinations of generation indices on A, 
there should be much stronger bounds because leptoquarks could mediate interactions 
forbidden in the Standard Model by lepton family number conservation and the absence 
of flavour-changing neutral currents (FCNC). 

In this paper, we concentrate on bounds from existing accelerators, meson-anti-meson 
(KK, DD, BB) mixing, rare lepton and meson decays, and a few electroweak tests. Most 
of these, have, of course, been previously calculated in various models, |5], [7], |TI| [TTJ, 



12] , g, [L3|, [L4|, [L5|, |j so we wish to list the bounds in as complete and model independent 
a fashion as possible. We neglect bounds from CP violation, because these constrain the 
imaginary part of the coupling, not its magnitude. 

In the remainder of this section, we introduce the interactions we are allowing the 
leptoquarks to have. There is then a short review of some theoretical models in which 
leptoquarks appear, followed by sections on pre-HERA accelerator mass bounds, con- 
straints from meson decays, meson anti-meson mixing, muon physics and r decays, and 
a section of bounds that are none of the above. 

We will assume that the couplings are real, so we neglect any constraints arising from 
leptoquark contributions to CP violation. (This simplifying assumption was not made in 
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|7|, [13], [15], [16], [L7|, |IBfl.) We will discuss each calculation in the text, and usually list the 
bounds (without generation indices). Most of the numerical upper bounds are collected 
in the tables at the end. 

We neglect QCD corrections to all the rates that we calculate, except b — > 57. This 
may be a poor approximation, but we would not expect this to change our bounds on 
A 2 by more than a factor of 2. 

There are seven renormalizable B and L conserving quark-lepton-boson couplings 
consistent with the SU(3)xSU(2)xU(l) symmetries of the Standard Model for both 
scalar and vector leptoquarks, and each coupling carries generation indices for the two 
fermions (these are suppressed in equations ([I]) and (fj)). The scalar and vector interac- 
tion lagrangians are therefore 



and 



C-s = { ^LS o q c L iT 2 £L + ^RS u c R e R )Sl + \ RS ~ o d c R e R Sl + (\ L s 1/2 Ur£ l + 

^RS 1/2 qLir 2 e R )S\ /2 + ^ L ~ Sl/2 d R £ L S\ /2 + X LSl q c L ir 2 r£ L ■ S\} + h.c. 

Cv = { (Xlv c1l1/l + ^Rv dRln e R) V o ] + XrvUrI^rV^ + (\ L v 1/2 d%rr/L 
+XRV 1/2 q C Lln e R) V i/2 + X lv 1/2 u c r1/lV{ 1 / 2 + \lvAl1^£ l ■ Vf } + h.c. 



(2) 



where the SU(2) singlet, doublet and triplet leptoquarks respectively have subscipts 0, 
1/2 and 1, and the L/R index on the coupling reflects the lepton chirality. We have 
written these lagrangians in the so-called "Aachen notation". In the first column of 
tables 3 and 4, we have written the interactions in this notation, and also {in curly 
brackets}, in the notation of reference 0. (In this alternative formulation, fermion 
number two leptoquarks have couplings g, fermion number leptoquarks have couplings 
h, and SU(2) singlet, doublet and triplet leptoquarks respectively have subscipts 1, 2 
and 3.) The fermions in our notation are chiral: here and in the tables, f R = (Rfyj 

and /;: (/,/•)•. 

Generation indices are superscripts and the lepton family index comes first: A 2 - 7 cou- 
ples a leptoquark to an ith generation lepton and a jth generation quark. If leptoquarks 
also carry generation indices, cancellations between the different generations of lepto- 
quark are conceivable (vaguely similiar to the GIM mechanism). In this case our bounds 
would apply to the sum over leptoquark-induced amplitudes, but not to the coupling 
constants and masses of each generation of leptoquark. We ignore this possibility. In the 
sections on rare meson decays, we quote bounds on e 2 (see (|T5| ) and (|T6|)), defined as A 2 
for vectors and A 2 /2 for scalars. 



2 Theory 

We present here a short review of some extensions of the Standard Model that could 
contain leptoquarks. The first possibility is Grand Unified Theories (GUTs) (19|], where 
leptons and quarks usually appear in the same multiplet. In consequence it is quite 
normal to have gauge and Higgs bosons coupling to a quark and a lepton (with or 
without B and L violation). Secondly, there are extended technicolour theories |[24|, 
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where quarks and leptons individually appear in multiplets of the dynamically broken 
extended technicolour group. The other particles in each multiplet are new fermions, 
that would appear at low energies in fermion-anti-fermion bound states, some of which 
are leptoquarks. And finally in substructure models p5| , the "preons" in a quark and 



lepton could combine to form a scalar or vector leptoquark. 



2.1 Grand Unified Theories 



Leptoquarks first appeared in Pati and Salam's SU(4) model p0|, where lepton number 
was treated as a fourth colour: the four weak doublets of each generation are arranged 
as a four of SU(4). SU(4) is then spontaneously broken, so that the gluons remain 
massless and the leptoquark gauge bosons become heavy. In this model, the leptoquarks 
induce flavour-changing neutral currents and lepton family number violation; they can 
also violate B, but their contribution to proton decay is severely suppressed ||2 1|| . 

Constraints on Pati-Salam type leptoquarks have been discussed in [|0j and |7|, ||. 
One would expect these vector leptoquarks to have full-strength couplings to a lepton 
and a quark of the same generation, giving, for instance, a large (~ X 2 /mf q ,X ~ 1) 
contribution to K — > efj, (s + \i — > d + e) . They could also interact with quarks and 
leptons of different generations, but one would expect this to be Cabbibo suppressed. So 
the best bounds on these leptoquarks should come from interactions involving lepton- 
quark pairs of the same generation. SU(5) [E2J GUT models contain vector leptoquarks, 
but these acquire grand unification scale masses when the gauge group breaks to SU (3) x 
SU (2) x i7(l), so are of little interest to us. This is a generic feature of GUTs: the quarks 
and leptons share multiplets, so there are necessarily leptoquark gauge bosons; however 
the lepton-quark symmetry is broken at the GUT scale, so these gauge bosons have large 
masses. (Note, however, that these are not neccessarily of order 10 15 GeV; the "GUT- 
scale" for Pati-Salam could be ~ 10 5 GeV.) SU(5) also contains scalar leptoquarks in 
the Higgs sector. The Standard Model Higgs doublet that breaks SU(2) shares an SU(5) 
5 with an SU(2) singlet leptoquark. Given the fermion representation in SU(5), this 
leptoquark necessarily has 5-violating interactions because it couples to two quarks 
and also to a lepton and a quark. The scalar potential must therefore be arranged to 
give these leptoquarks very large masses, in which case they are of no interest to us. 
Alternatively, one can break SU(5) invariance by setting the leptoquark-quark-quark 
couplings to zero, as was done in reference |ll|], which allows all five components of the 



scalar 5 to be light. It is of course possible to have scalar leptoquarks in SU(5) with B 
and L conserving interactions by putting them in a representation other than the 5. An 
example of this can be found in J^, where it is shown that the minimal SU(5) model plus 
a light (mi q ~ 100 GeV) SU(2) doublet leptoquark would be compatible with present 
experimental results on proton decay and sin 2 #w This B and L conserving leptoquark 
(= 5*1/2, see table 3) occupies the same position in the SU(5) 10 as the quark doublet, so 
its interactions with light Standard Model fermions are of the form cIrIlSI/2- The other 
members of the 10 are assumed to be very heavy. 

The field theory limit of the heterotic superstring is some "GUT-like" gauge field 
theory, which is broken at or below the compactification scale to the Standard Model 
and possibly some extra {7(l)s: SU(3) x SU(2) x {7(l) n+1 . This low energy limit may 
contain leptoquarks for the same reasons that GUTs do. For instance, certain Calabi- 
Yau compactifications have G = E e /H as the four- dimensional gauge group below the 
compactification scale, where H is a discrete symmetry group that can be chosen such 
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that G is the Standard Model (xU(l) n ). The low energy fields are then expected to sit 
in the 27 of E 6 , 

and for each generation consist of quarks and leptons, a right-handed neutrino, the 
two Higgs doublets of the Minimal Supersymmetric Standard Model, an extra SO(10) 
singlet, and a pair of SU(2) singlets, which may have thr couplings of either diquarks or 
leptoquarks, but not both simultaneously ||13|| , as this would destabilize the proton. The 
constraints on leptoquarks in such superstring-derived models were calculated in 
and are extensively reviewed in 



II, 



rj. 



2.2 Technicolour 

If our knowledge of quantum field theory is extrapolated to energies far above those 
experimentally accessible today, one discovers that the masses of fundamental scalars 
are generically the same size as the cutoff energy. So to keep the mass of the Higgs at 
the electroweak scale, one must either fine-tune the parameters of the theory to many 
decimal places, or introduce new physics that explains the comparatively small scalar 
mass. Both supersymmetry and technicolour |S| do this. 

In a technicolour theory, all scalars are bound states of fermion anti-fermion pairs, 
like the mesons of QCD. One introduces new electroweak doublets and singlets (tech- 
nifermions) which are multiplets of a non-abelian gauge interaction called technicolour. 
The new interaction becomes strong at the electroweak scale, and, by analogy with QCD, 
one expects a techni-fermion condensate to form, breaking some of the global symme- 
tries of the theory. Three of the goldstone bosons arising from this global symmetry 
breaking become the longitudinal components of the electroweak gauge bosons. This 
neatly gives masses to the W and the Z, without a fundamental scalar Higgs. In fact, 
if one temporarily ignores the problem of generating fermion masses, one can simply 
introduce an electroweak doublet of left-handed technifermions, and their right-handed 
singlet counterparts. This has an SU{2) L x SU{2) R global symmetry if one neglects the 
(comparatively weak) standard model interactions. When the technifermion condensate 
forms, the global symmetry is assumed to be broken to SU(2)y, and the three goldstone 
bosons associated with the broken axial currents (~ pions in QCD) are eaten by the 
electroweak gauge bosons. 

To be a viable substitute for the scalar Higgs, technicolour must also produce fermion 
masses; this is done in extended technicolour (ETC) models. The idea is to introduce 
more technifermions, and gauge interactions involving ordinary and technifermions that 
are broken at very high energies. The gauge bosons of these new interactions (which we 
refer to as ETC gauge bosons) presumably acquire mass by some dynamical symmetry 
breaking mechanism of their own, and afterwards induce four-fermion interactions in- 
volving two technifermions and two ordinary fermions. After the formation of the tech- 
nifermion condensate, the four-fermion interaction becomes a "standard model" mass 
term. 

It is desirable to introduce a whole standard model generation of technifermions so 
that the ETC gauge bosons do not have to carry SU (3)xSU (2) xU(l) quantum numbers. 
This means that there will be "leptoquark technimesons" made up of a techniquark 
and an anti-technilepton. It is easiest to see why leptoquarks naturally arise in ETC 
models by looking at a simple example. An extra generation, plus right-handed neutrino, 
who each transform as an N under the technicolour SU(N), are added to the Standard 
Model, and one removes the Higgs sector. One also must add other interactions and 
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particles (including the ETC gauge bosons) at higher energies. The new technifermions 
are assumed to interact with each other via the techniforce, and with other Standard 
Model particles via colour and electroweak interactions, and through the effective four- 
fermion vertices induced by ETC gauge bosons. Technicolour becomes strong around 1 
TeV, a technifermion condensate is assumed to form, and some of the global symmetries 
of the technicolour sector are dynamically broken. 

In QCD, if one neglects the masses of the u and d quarks, there is a global SU{2)i J x 
SU{2) R symmetry. This is broken by the condensate to SU(2) V , and the pions are the 
(pseudo-)goldstone bosons of the broken axial symmetry. Similar behaviour is expected 
in a technicolour theory; in the model considered here, if one neglects electroweak, colour 
and four-fermion interactions, there is a global SU(8)l x SU(8)r symmetry (assuming 
complex representations under the technicolour group) which is expected to break to 
SU(8)v- This produces 63 goldstone bosons. Three will be eaten by the electroweak 
gauge bosons, and the others acquire masses due to the neglected interactions. It is 
clear that among the 60 goldstone bosons, there will be "leptoquark mesons" composed 
of a techniquark and an anti-technilepton, bound together by the technicolour force. 
This scalar leptoquark will interact with a quark and an anti-lepton carrying the same 
Standard Model quantum numbers (it would be easy to construct a model where the 
leptoquark had fermion number 2). 

The coupling of a technimeson leptoquark to a lepton and a quark occurs when the 
technifermion constituents of the leptoquark exchange an ETC gauge boson and turn into 
an ordinary quark and lepton. The "Yukawa" coupling can only be predicted in a specific 
model, but it is expected to be of order mj/A or m q /A, where A is the technicolour scale. 
Constraints on leptoquarks in technicolour models were computed in [0, §, []]. 



2.3 Composite Models 



suns! 



Leptoquarks frequently arise in composite, or substructure, models [go 
In these scenarios, the leptons and quarks, and sometimes the gauge bosons, are 
assumed to be bound states of more fundamental particles. In principle, one could 



hope that this would explain the quark and lepton spectrum, and reduce the number of 
fundamental fields (as quarks did for the hadronic spectrum). However, at the moment, 
most models do not require fewer constituents ( "preons" ) than there are particles in the 
Standard Model. 

One of the fundamental questions in building theories where the Standard Model 
particles are composite, is to explain why their masses are much smaller than the com- 
positeness scale. A solution for the fermion masses is to have them forbidden by a 
global chiral symmetry in the strongly interacting sector of the theory. This symmetry 
is broken by the Yukawa and SU(3) c x U(l) couplings, so one would expect the masses 
generated by these symmetry-breaking effects to be small. Models depending on preon 
chiral symmetries to preserve the masslessness of composite fermion bound states must, 
however, satisfy the 't Hooft anomaly matching conditions ||30|| . These require that 
chiral anomalies calculated in the effective composite theory match those calculated in 
the preon subtheory 



30, 31, 32 



A solution of the anomaly matching conditions is to 
reinterpret the "tumbling" gauge theories of dynamical symmetry breaking as com- 
posite models satisfying the 't Hooft conditions via the principle of complementarity. 
Although solutions of the anomaly matching conditions have been pursued by a variety 
of methods, (for a review, see p5f), the construction of a fully realistic model remains 
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an open problem. Another method for keeping the composite fermions light is to make 
them the supersymmetric partners of goldstone bosons in a theory with a spontaneously 
broken approximate global symmetry. Their masses must then remain small even af- 
ter the breaking of supersymmetry and the global symmetry. It is even more difficult 
to explain how composite gauge bosons could have masses much below the substructure 
scale; the Case-Gasiorowicz- Weinberg- Witten theorem forbids the presence of interacting 
massless spin-1 composite particles, suggesting that all composite vector bosons should 
have masses of order the compositeness scale |36] . 



If the quarks and leptons are composite particles, they should have form factors. The 
"substructure contribution" to g — 2 of the muon is expected to be of order (m^/Ac) 2 , 
where A c is the compositeness scale. This constrains A c <; 1 TeV, which makes it difficult 
to associate the weak scale with the substructure scale fas was done in 1211) in realistic 



models. There are a number of stronger, but more model-dependant bounds on A c |]38 

Leptoquarks appear naturally in many substructure models|^[. This is not sur- 
prising: if a constituent particle of a Standard Model fermion carries quark or lepton 
number, then a composite quark could turn into a lepton in the presence of a compos- 
ite lepton (that turned into a quark) by exchanging the appropriate constituents. The 
bound state consisting of the exchanged constituents would be a leptoquark. This sug- 
gests that leptoquarks in composite models could naturally induce interactions where 
flavour or generation number was conserved overall, but separately violated in the quark 
and lepton sector. 



3 Mass Bounds 

There have been numerous searches for leptoquarks at existing accelerators. Possible 
searches at ee colliders were discussed by Hewett and Rizzo p7| who point out that the 



measurement of the R-ratio from PEP and PETRA constrains all scalar leptoquarks to 
have mi q > 15 GeV. The leptoquarks are assumed to be pair-produced in the decay of 
a virtual photon, so this bound depends only on their electric charge which is known to 



be > 1/3 (see tables 1 and 2). The same measurement from AMY |4(| gives 

mi q > 22.6 GeV (for scalars). (3) 

We expect a similiar bound for vectors 

mi q > 20 GeV (for vectors). (4) 

These are the best model independant lower bounds on the mass (that we are aware 
of). There have also been leptoquark searches at LEP and the pp colliders, but these 
look for leptoquark decay products, so only apply if the leptoquark-lepton-quark Yukawa 
coupling is "sufficiently" large (although this is not a significant constraint; see below). 
The LEP lower bound |^T| is 

m lq > 44 GeV (5) 

for scalars, and a similiar bound should apply to the vectors. This was calculated assum- 
ing that the Z Q turned into a pair of leptoquarks, which then decayed to jets and two 
leptons (ee, rf , vv; the mass bound actually depends slightly on the type of lepton 
f4l| — (|) is an approximation). Since all the leptoquarks are charged, they all couple 



to the Z D , although this is rather weak in the case of the Q em = 1/3 scalar singlet (see 
tables 1 and 2). 
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Delphi has also looked for Z decays into an on-shell and an off-shell scalar lep- 
toquark. They conclude that 

m lq > 77 GeV (A > e) (6) 

for leptoquarks decaying to a quark and an electron or a muon. Note that unlike the 
bounds (H) and (|9]), this only applies for large (A > e) leptoquark-lepton-quark couplings. 
In principle, leptoquarks could contribute at one loop to the decays Z Q — * £i£ 2 , with 
l\ 7^ £ 2 , but the rates for such processes are so small that LEP would not expect to see 
them even if the leptoquark yukawas were of order 1 |42| . 



The decay rate of a leptoquark to a lepton and a quark is 

T = lq (for scalars) (7) 

l07T 

r = ^rj-^ (for vectors) (8) 

247T 

so one would need the coupling constant to be very small (^ 1CT 8 ) for the leptoquark 
to decay outside a detector. Such a small renormalizable coupling would be unexpected, 
so it should be safe to assume ra iq > 44 GeV. (Note that equations (^) and (§) give the 
leptoquark decay rate to a specific flavour of quark and lepton. The total decay rate 
would be the sum over all types of kinematically accessible quarks and leptons.) 



The UA2 bound |44| on scalar leptoquarks decaying to e+ jets is mi q > 67 GeV, and 
CDF has a bound of 113 GeV (95%) f45f . The strongest pp collider bound is a DO result 
| 46fl for leptoquark pairs decaying to two electrons plus jets: 



mig > 126 GeV (lq -> eg) (9) 
where the leptoquark is assumed to decay entirely to electron plus jets. The cross section 



for vector leptoquark production is being calculated |47|], so bounds on vector leptoquarks 
from CDF and DO should be possible in the future. 

It is not clear that leptoquarks that could be seen at HERA must satisfy these bounds; 
although they must be produced in a quark-electron (or positron) collision, they could 
have stronger couplings to muons or taus, decay to fi or r plus jets, and thereby avoid 
this constraint. We will nonetheless use mi q = 100 GeV in our numerical results. 

The ZEUS and HI Collaborations at HERA have (mass-dependent) limits on the 
couplings of scalar and vector leptoquarks produced in an ep collision and decaying to 



e + X or v + X [48]. The exact bounds depend somewhat on the type of leptoquark, 



and will of course improve with time. The present limits rule out roughly from 

A > .05 m lq = 25 GeV (10) 

to 

A > 1 mig = 225 GeV (11) 
for the SU(2) singlet and triplet leptoquarks (scalar and vector), and from 

A > .1 m lq = 25 GeV (12) 

to 

A > 1 m lq = 150 GeV (13) 
for the doublets. Note that these are just very rough approximations to the present data. 



See |48[ for accurate plots. 



8 



4 Meson Decays 



Since leptoquarks have not been directly observed at any high energy colliders, one 
must look for traces of them in other interactions (meson decays, for instance). At 
energies well below their mass, the leptoquarks can be "integrated out" , leaving effective 
four-fermion vertices involving two leptons and two quarks, and a photon-lepton-lepton 
(or quark-quark) magnetic moment vertex. We will come back to the latter in sections 
6 and 7. The four fermion vertices contract the spinor indices of a lepton and a quark, 
but can be Fierz-rearranged to contract the quark indices together; in this form one 
can estimate the quark matrix elements between initial and final meson states. The 
leptoquark induced effective four-fermion vertices are listed in tables 3 and 4. The 
hermitian conjugates of the listed operators are of course also possible. Most of the 
leptoquarks couple to quarks of only one chirality, so the Fierz-rearranged vertices must 
be of the form V ± A, similiar to the weak interactions. The only exceptions to this are 
S , S1/2, Vq and Vy 2 , which have two interactions: a "left-handed" one with coupling 
constant Xl, involving the leptoquark, a left-handed lepton and a quark, and a "right- 
handed" interaction coupling the leptoquark to a right-handed lepton and a quark. The 
Fierz-rearranged effective four-fermion vertex arising when these two interactions appear 
at opposite ends of the leptoquark propagator contains tensor, scalar and pseudoscalar 
pieces. We will ignore the tensor operators, because their matrix elements are difficult 
to estimate, and just list the scalar ± pseudoscalar parts. Some experimental results 
will give much better constraints on the coefficients of the S ± P matrix elements than 
the V ± A ones. In these cases, we will list bounds on the products Xl^r. In general, 
however, we will avoid the "mixed coupling" constraints because they do not apply if 
one of the interactions does not exist. 

In the Standard Model, many mesons are forbidden from decaying via the strong or 
electromagnetic interactions. Since such mesons decay at a rate consistent with the weak 
interactions, one would roughly expect the effective dimensionful coupling ~ \ 2 /mf q of 
the leptoquark four-fermion vertices to be as small as that of the weak interactions: 

where V is the appropriate CKM matrix element. 

However, we can get much stronger constraints on the leptoquark couplings that 
induce interactions not present in the Standard Model. A variety of meson decay con- 
straints were considered in references |fj, [| [7], [| [IT], O, [T3|, |T3|. The first "non- 
Standard" type of leptoquark-induced behaviour we will consider is the lack of chiral 
suppression in meson decays due to pseudoscalar matrix elements. The V-A vertices 
of the weak interactions give a matrix element squared for the leptonic decay of pseu- 
doscalar mesons proportional to EiE n — k 2 (where Ei,E n and k are the energies and 
momentum of the outgoing leptons), because the helicity (~ chirality for a relativistic 
particle) of one of the outgoing leptons must be flipped via a mass term to conserve 
angular momentum. The decay rate is therefore strongly suppressed for relativistic lep- 
tons. There is no such chiral suppression in meson decays induced by leptoquarks with 
both left and right-handed couplings (pseudoscalar matrix elements, see below). 

One would also expect leptoquarks, if they exist, to induce lepton family number 
violating decays; even if the leptoquark Yukawa interactions are "generation diagonal" 
in the sense that they only couple first generation quarks to first generation leptons, this 
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would still allow, for instance a kaon to decay to e\i. It would similiarly be very natural 
for leptoquarks to induce flavour- changing decays suppressed in the Standard Model. 
Leptoquarks can induce effective four fermion vertices of the form 

tj TfL'fl 

^^(gV^VXrv'O (is) 



and/or 



t] van 



e J e 



™? lq 



-(q j P q q n )(£ m P e t) (16) 



(see tables 3 and 4) where P q and P are chiral projection operators = L or R, and mi q 
is the leptoquark mass. The coefficients of the effective operators in the third column 
of tables 3 and 4 are generally twice as large for vectors as for scalars. In equations 
(rj) and ( IBD we have introduced a coupling e 2 = A 2 (e 2 = A 2 /2) for vectors (scalars); 
the constraints we compute on e 2 will therefore be bounds on the square of the vector 
coupling, or half the square of the scalar coupling. As is clear from tables 3 and 4,V±A 
quark currents are induced by leptoquarks with couplings of the same chirality at both 
ends of the propagator. The effective coupling of V ± A 4-fermion vertices is therefore 
e 2 L /m 2 q or e 2 R / m 2 q . In all the following, we will quote bounds on e 2 L with the understanding 
that the constraint also applies to e 2 R . (Remember that the coupling constant subscript 
applies to the lepton chirality.) The scalar and pseudoscalar vertices are induced by the 
leptoquarks S a , Si/ 2 , and Vy 2 with couplings of opposite chiralities at either end of 
the propagator. The effective coupling of the S ± P vertices is therefore e^eR/m^. 

The matrix element squared for the decay of a pseudoscalar meson M (qiq n bound 
state) to the leptons t and I m is 

\M\ 2 = Mi(4er) 2 + (4eT) 2 }<iA 2 (E l E m -k 2 ) 

lq 

+ MeT? + {eleTf\P\^E m + k 2 ) (17) 
+2(4 ef n + 4e R m )PAm M {4eTE i m m - 4e r E n E i m m )} 

where 

P = - < >~ ——————— mM (18) 

2 2 m q j + m q n 

and 

Ap^ = ^ < 0|gV7V|M>= (19) 

(p M is the meson 4-momentum and J'm the decay constant). The approximation for P is 
a current algebra result, so not appropriate for the heavier mesons; we will use it anyway 
to get qualitative results. Equation(|TTD gives a meson decay rate 

r(M-« m ) = ^Mi^erf + ieleTfWMA^E^-k 2 ) 

M lq 

+ [(4eT) 2 + (4eT) 2 }P 2 (E*E m + k 2 ) (20) 

+2(4 er + eleT)PAm M (e l leTE l m m - E imm )} 

where k is the magnitude of the lepton 3-momentum in the centre-of-mass frame, and 
the absence of chiral suppression in the pseudoscalar matrix element is clear. 
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4.1 Pion Decays 

In the Standard Model, charged pions decay principally to jiv because the decay to ev 
is suppressed by angular momentum conservation, as discussed in the previous section. 
The experimental ratio 



R exp = = 1.231 ± .006 x 10" 4 (21) 

r(7T+ — > fiv) 



agrees with the Standard Model prediction |50| 



Rth = K rn 2 e f + L235 ± 0Q4 x 1Q _ 4 
(m£ — wit) 

where A is a radiative correction. Using ( PH| ) to compute the contribution to this ratio 
from the interference of the Standard Model axial vector amplitude and the presumably 
small leptoquark-induced effective pseudoscalar operators, one can require 

P l e 2 e R m ?r e L e fe m n ^ \Rexp — Rth\ 



m n A V2G F mf q m e y/2G F mf q m^ R th 

If we assume that there is no spurious cancellation between the pseudoscalar contri- 
butions to 7r — > fiv and 7r — ► ev (first and second terms of equation (|23"D ), then this 
gives 



where we have used the current algebra result P ~ 7 ' j^m^ @, and n is a light neutrino 
index (m^ <C m e ). 

We get bounds on (e^/m^) 2 from the interference between W and leptoquark ex- 
change, which, if small, must satisfy 

where n is again an arbitrary light {m v <C m e ) neutrino index. Assuming that leptoquark 
couplings to muons and electrons are not equal (e^ 1 ^ e^ 1 ), we can conclude that: 

«-«< 2 » r '(i»^) ! • < 26 > 

If the leptoquark coupling is independant of the lepton generation, so that e 11 = e 21 , 
then leptoquarks will not affect the R ratio. However, (e 11 ) 2 [= (e 21 ) 2 in this case], is still 
constrained by quark- lepton universality (see section 8.1), which gives a better bound 
than this one. 

Equations fl24]) and (|2~6"D give bounds on the couplings of the effective vertices of tables 
3 and 4 that involve u, d, e (or //) and any flavour of light (m u <C m e ) neutrino. We list 
the exact constraints in the tables. Since the notation in the tables is rather obscure, we 
will discuss in detail the first four rows of the table listing bounds on vector leptoquark 
couplings (table 5). The four-fermion vertices have an effective coupling constant e 2 /mf q , 
so the bounds on e 2 depend on mf q . We take mi q = 100 GeV to compute the numbers 
in the tables, so the numerical bounds scale as (m^/100 GeV) 2 . The rows are labelled 
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by experiments, and the columns by coupling constants; an entry in the table is the 
bound on a particular coupling constant from the given experiment. These constraints 
are in general generation dependent, so we list the coupling constant generation indices 
in parentheses. For example, we have just found that the upper bound on R constrains 
axial-vector quark operator couplings to be less than 2 x 10 _3 (m; g /100 GeV) 2 . Only 
leptoquarks inducing an interaction between u, d, v and e or ji could contribute to R, so 
we see from table 4 that we have bounds on X 2 LVo and A| yi . Since the neutrino in the 
decay n + — > ev could be any species providing that it is light, we have the bounds 



AiV \f v < 2 x 1(T 3 ( m " ) (27) 
LV ° LV ° V100 GeV/ V ; 

*lv, < 2 x 10- 3 ( mi « ) 2 (28) 
LVl LVl V100 GeV/ V ; 

where n is the generation index of a light neutrino. This is the first row of the table. 
From (|26|), the same bounds apply to leptoquark couplings inducing the decay 7r — > /iz/, 



which gives the second row of the table. We also have bounds on pseudoscalar operators 
from the R ratio. Again looking in table 4, we see that there are such operators with 
effective couplings ^LVo^RVo/ 171 ^ and \ LVl/2 \ RVl/2 /mj q . Checking carefully in the second 
column to get the right generation indices on the right coupling constants, we see that 
the bound ( f2~4"D gives 

« < 5 x 10- (29) 

^^ 1/I <5xlO-( I i) 2 (30) 

as listed in the third row of the table. Note that bounds that cross two columns apply 
to the product of the couplings in those columns, with the first set of generation indices 
applying to the first column coupling constant and the second set to the second column. 
The fourth row of the table is the constraint from (^4j) applied to pseudoscalar operators 
involving muons. 

From the upper bound BR(7r + — > Jiv e ) < 8 x 10~ 3 , one can constrain leptoquark 
induced effective interactions of the form 

2 

^{u^Pd){e llx Lp) , e -^(uPd){eLv) (31) 

m lq m lq 

to satisfy 

^<.25GV, ^<5xlO- 3 G F ^ (32) 
mf q mf q m n 

where we have again used P = 7 ' j^m^. 

There is also an upper bound on the lepton flavour violating decay n° — > / u ± e =F : 
BR(7r° — > /U ± e =F ) < 1.6 x 10~ 8 , and the branching ratio for tt° —* e ± e =F has been measured 
f51"| to be ~ 7 x 10 -8 . Unfortunately, since the neutral pion decays electromagnet ically 
7r° — > 77 with a lifetime eight orders of magnitude smaller than that of the charged 
pions, these do not translate into very strong bounds on the leptoquark couplings. We 
will get much better constraints, with the same generation indices, from muon conversion 
on nuclei, and atomic parity violation. 
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4.2 Kaon Decays 

Charged kaons decay to fiv^ (BR = 63.5 %) and a variety of other final states, most 
involving pions. The decay to ev e is suppressed by angular momentum conservation as 
in the pion case, so if we require that the leptoquark contribution to the branching ratio 
be less than the experimental error, we get the results listed in the tables. We assume, 
as in the pion case, that there is no cancellation between the leptoquark contributions 
to K fiv and K — > ev, and we use P ~ to estimate the pseudoscalar 

contribution. 

If kaons decayed via a leptoquark, it would be just as natural to have \w e in the final 
state as fiv^. The upper bound on the ratio |38| 



BR(K -> pup) 
implies 

el (e L e R ) 2 P 2 
%mj q G 2 F \V su \ 2 ' 2m\ q G\\V su \*p K ml 



^-^«)< 6 .3xl0-3 (33) 



< 6.3 x 10~ d (34) 



where we again use P ~ mx/if/2 to calculate the numerical constraints listed in the 
tables. 

There are also bounds from the absence of the flavour-changing decays K — > n£i£2, 
where £i£ 2 = ee, \ip, fie or emu. The Standard Model does not allow a strange quark to 
turn into a down quark at tree level (FCNC), whereas leptoquarks could easily induce 
an sd££ vertex. If we assume, using isospin symmetry, that 

< K + \0\n° >= -L < K + \6\ti + > (35) 
V 2 

(where O is some isospin 1/2 operator) and neglect all lepton masses, we can require 

4 < iGM BR(K^«J) (36) 

giving the constraints listed in the tables. Neglecting the masses should introduce an 
error of less than a factor of 2, because BR(K — > n ev) ~ 1.4 BR(K — > -K°pv). 

Finally we consider constraints from leptonic decays of neutral kaons. Most scalar and 
vector leptoquarks would allow the decays Kl — * pp, ee, fie, which are suppressed in the 
Standard Model by the absence of flavour-changing neutral currents (and lepton family 
number violation). There are both axial vector and pseudoscalar quark matrix elements 
that can contribute to these decays. We will neglect the pseudoscalar contributions in 
the decays to fie and \ip because they constrain the product eien, and the bounds on 
the individual coupling constants {e 2 L ,e 2 B ) are only a factor of m^/m^- weaker. Using 
equation (|20| ) with A = f K /2 we find that 

T(K L - eV) - (m *~? ;) ' 4/X < L2 x lir " GeV < 37 > 



QAnm 6 K m\ 



which implies 



e 2 T < 6 x 1(T 7 
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One can similiarly compute the decays Kl — > fJ>fi and Kl — > ee from equation ( p0|) . This 
gives constraints as listed in the tables. We have included the bounds from effective 
pseudoscalar operators in the decay to ee, using P ~ because they are consid- 

erably stronger (not suppressed by the small electron mass). Note, however, that scalar 
leptoquarks do not induce pseudoscalar operators of the form (d> Pd rn ){e l Pe m ) (see table 
3), so these bounds only apply to Xlv ^rv an d Xlv 1/2 ^rv 1/2 - 



4.3 D decays 

Since the mass of the D + meson is much greater than that of the muon, the decay 
D + — > /i + z/ is suppressed by angular momentum conservation, like tt + — > e + u. The 
present upper bound on the rate for D + — > fi + u is roughly a factor of 2 larger than 
the Standard Model prediction, so is not the best place to get bounds on leptoquarks 
that induce axial vector effective vertices. However, it does strongly constrain effective 
pseudoscalar quark vertices involving c, d, \x and any neutrino: 



{m 2 D -mp 2 p 2 (e L e R ) 



^ P 2 ^p- < 4.5 x 1(T 16 GeV (39) 



or, (assuming P ~ f D m D /2, with |52| f D ~ .2 GeV) 

«*> < 6 x irS imkv) 2 ■ (40) 

A much weaker bound can be calculated for the effective coupling of vertices involving 
c and s quarks, a muon and a neutrino, from the upper bound BR(D$ — ► /i + z/) < .03. 
We do not list this in the tables. 

The CKM matrix element V c d is measured in AC = 1 neutrino interactions (y + d — > 
fi + c) to be .204 ±.017 | 38fl . Requiring that leptoquarks contribute less than the total 



observed rate (a conservative assumption) gives 

el < ^\V cd \m 2 q (41) 

for vertices involving c, d, \x and v^. The upper bound on the branching ratio for D + — > 
[i + v gives slightly weaker constraints, but applies to all flavours of neutrino, so is included 
in the tables. 

The CKM matrix element V cs is determined to be 1.0 ± .2 from the decay D° — > 
K~e + u, so one can require 

e? 4:Gf , ln . 

where e\/mf q is the effective coupling for an axial vector vertex involving c, s, e, and any 
neutrino. 

We can constrain the effective coupling of the flavour- changing operator {u^^Pc) (l^Pl) 
(I is an electron or a muon) from the upper bound on the branching ratio BR(D + — > 
7T + U), if we assume, by isospin symmetry, that 

< D + \u^Pc\n + > = < D°\d'fPc\'K- > (43) 



14 



(where P = L,R). The D° is observed to decay to tt ev with a branching ratio of 3.9 
xlCT 3 , so if we neglect all lepton masses we can require 



p 2 

L 



V2 



mf q 4:G F \V cd \ 



< 



r D o BR(D+ -> Ti+ll) 
\| r D+ BR(D° -»• ix-ve) 



(44) 



which gives the constraints listed in the tables. 

As in the case of kaons, one expects leptoquarks to induce flavour- changing D° decays 
to lepton pairs (pfi,ee,pe). Most of the leptoquark interactions induce effective axial 
vector quark and lepton currents, and since m D 3> m 



i-i- 



m P 



the contribution of these 
currents to the decays D — > pfj,, ee, pe is suppressed by angular momentum conservation. 
One gets strong bounds on the coupling constants of effective pseudoscalar operators 
from the leptonic decays, but better bounds on the axial vector operators from the 
previously considered D + — > decays. Using (p0|) to calculate the rates for D° decay 
with the approximation flUf), the limit BR(D° — > fie) < 10~ 4 [SS| gives 



e L e R < 4 x 10" 



miq 



100 GeV 



(for D° -> fie) 



(45) 



Similiarly for BR(D° -»• up) < 1.1 x 10" 5 |5S[ one gets 

2 



e L e\ < W 



mi q 



100 GeV 



(for D° -> (ip) 



(46) 



and from BR(D° -> ee) < 1.3 x 10" 4 |8 



e L e^ < 4 x 10~ 3 



100 GeV 



(for £>° -> ee) 



(47) 



4.4 B decays 

One can get constraints on the leptoquark couplings to b quarks from the upper bounds 
on the flavour- changing decays B — > ttX (££ — fj,p,ee), B — > £i! 2 K = fj,p, fu,e,ee), 
and B° — > pfi, ee, pe. If we approximate the hadronic matrix element in the B t\t<iK 
decays via heavy quark effective theory |54||, we get some of our best constraints from 
these processes. We nonetheless calculate and list the bounds from B — > flX and 
B° — > first because they should be more dependable. 

The upper bounds on the leptonic branching ratios are stronger than those for the 
inclusive decays B — > ££X, but the constraints on the couplings are weaker because 
the theoretical rates are suppressed by angular momentum conservation (small lepton 
masses). We therefore list, for the leptonic decays, bounds on the coupling constants 
of effective pseudoscalar vertices, and bounds on coupling constants whose generation 
indices would allow B° — > pe, because the latter are not constrained by the inclusive 
decay data. 

The decay B~ — ► evX (B~ — > fiuX) is observed to have a branching ratio of 12 % 
(11%). This is the rate expected in the Standard Model with \V c b\ = -05, so we can 
require that the leptoquarks mediating this interaction satisfy 



m 



V2 



cb\ 



(48) 
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The branching ratio for B — > rz/X has recently been measured to be 4.2 + . 72 — .68±.46% 
by Aleph []55| . One would expect the rate to be suppressed by about a factor of 2.5 with 
respect to B — > \wX by kinematics, so leptoquarks coupling a tau, a neutrino, a b and 
a c or u quark must also satisfy fl48|) . If we neglect differences in the kinematics between 
B + — > ££X and 5 + — >■ z/£X, where £ = e, [i, we can also constrain 

^ < 8 ^ |Kfe| 5i?(5+ ^ (49) 

(IjD and fl4~9|) give bounds on four fermion vertices involving ev b c (or eu b u, but we 
will get better bounds on this from (|50|), ££6s and iibd. 

Taking = .06, we obtain the constraints listed in the tables. The ratio of CKM 
matrix elements V u b/V c b is measured to be [TO .10 ±.03 by fitting the lepton energy 



spectrum in semi-leptonic B decays. We can therefore constrain leptoquarks coupling a 
b, u, e and a v to satisfy 

el , AG F 



m 



2 



&—i=r\Vvb\ • (50) 

Hq V2 

We can get constraints on four fermion interactions involving b, r, a down-type quark 
(d or s) and another lepton (I = e, /i or r) from the decays 5 — > £lX because the r 
decays very rapidly. Suppose first that £ = e, and the leptoquark mediates the decay 
B — > reX. We are interested in the r decaying to v + hadrons, in which case this 
would contribute to B — > ueX, or the tau decaying to euv, giving eevvX as the final 
state. The experimental limit on the branching ratio for B — > eeX is smaller than the 
measured rate for B — > veX, so we would get better constraints if we could assume 
B — > eevvX C B — * eeX. However, the CLEO Collaboration |57| required that the ee 
pair in its upper bound on B — > eeX come out back to back, and since the tau is not 
extremely relativistic, it is not clear that the electron from its decay would be at 180° 
from the e.We will therefore calcalate bounds from the rate for B — > veX. We assume 
that the phase space suppression due to the tau mass contributes about a factor of 1/2.5, 
so 

Tin v\ 1 e ' r(J? -> evX) 

- reX) = 8G|| ^ |2 . (51) 

We can therefore require that couplings involving b, r, e and an s or d satisfy 

^ < 2^,—-%—^ . (52) 



[Bi2(r -> z/X')] 1 / 



Hq 

If the 6 decayed to r, /i and a s or a d quark, and the r decayed to fivv, this would 
probably have been seen by UAl in their search |58| for B — > fip,X. We can therefore 
require 

±G F V cb I BR(B wX) 

mf q < V2 \ BR(B -> {ivX)BR(t -> ^ii/p) 1 J 

for couplings involving 6, r, /i and a or an s quark. 

Finally we consider the decays i? — > rrX. We assume that both taus decay to muons, 
and we include an approximation to the the phase space suppression due to the masses 
of the final state particles [} The effective coupling of four-fermion vertices involving a 

1 We thank Andrzej Czarnecki for creating a simple analytic approximation for us. 
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b, two rs, and a d or an s quark must therefore satisfy 



e 2 AGp 



m 



BR(B -> /i/iX) 



(54) 



where PS 1 is the phase space suppression factor: 



PS 



1 - 



2m T 



2.48 



(55) 



The rate for the decay B° — > eu, via an effective axial vector quark operator will be 
fl57|) with the obvious modifications: 



T(B° 



(m 2 B - ml) 2 



L 



f B ml < 3 x 1(T 18 GeV 



647rm 3 



B 



111 



(56) 



where we have used the recent CLEO bound BR(B° — > e/i) < 6 x 10~ 6 to get the 
numerical constraint, and we assume ~ .2 GeV to get the constraints in the tables. 
Using P ~ in ((20), we can estimate the decay via a pseudoscalar vertex to be 



.03^^ GeV 5 < 2 x 1(T 17 GeV 



m 



(57) 



which gives the bounds on Xl^r listed in the tables. These constraints, however, are 
only rough estimates. 

The most stringent experimental bounds on flavour- changing neutral currents involv- 
ing b quarks come from the exclusive decays B + — > X + £i£^ , where if £2 can be ee, fie 
or fip,, and X is a K or a 7r. However, to calculate a rate for this decay, we need to 
evaluate the matrix element < X + \s , ~f fl b\B + >. We do this by relating it at zero recoil 
to < X + \s , ~f fl c\D > via the heavy quark formulism and then making very simple 
approximations for the momentum transfer dependence of the form factors. 

The B to X matrix element can be written 



< X + \sYb\B + >- 



+ Px) 



f B - x (p 



B 



fx) 



(5? 



where f± are functions of the momentum transfer t = (ps — px) 2 - A similiar expression 
can be written for the D° — > X + matrix element, and the B and D form factors can then 



be related by the heavy quark symmetry 
and A /nib, this gives 

J -\- \ L max) J— \°max } 



Neglecting corrections of order m c /nib 




a(m c ) 



-1 -6/25 



f U ^ X it ) 
J 4. \ L max } 



(59) 



where t max is the maximum momentum transfer between the mesons {{tub — nix) 2 for 
the B decays, {m D — m x ) 2 for the Ds). 

The decay rate for B — > X£ + £~ is proportional to /? ||60|| , so if we were to ignore the 
t dependance of the form factors, and the mass of all the final state particles in both the 
B and D decays, we could require that 



m 



< SG 2 F \V\ 



lq 



2 m D m c 

m B nib 



u{m b ) 
a(m c ) 



12/25 



T{B^ 



XV 



T(D° -> Xlv) 



(60) 
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where |V| is the appropriate CKM matrix element, and T(B + — > X£ + £~) is the experi- 
mental upper bound on this rate. 

These approximations are rather extreme; the form factors can vary by up to an order 
of magnitude as t varies from i m j n to t max , and m^/mp is far from negligable. To model 
the momentum dependence of the form factors, we assume that they vary with t as 



/+(*) 



/+(0) 
1-t/m 2 



(61) 



(where we take the mass to be that of the decaying pseudoscalar meson for convenience; 
the vector meson mass might be more accurate, but for b and c mesons, these are very 
similiar), then ( |59D becomes (if we neglect mx/mB,mx/mi)) 



/*-*(0) 




a[m b ) 
a(m c ) 



-6/25 



/r x (o) 



(62) 



and the bound ([BOf ) is 



<8G%\V\ 2 ' 



m 



lq 



a{m b ) 
a(m c ) 



-i 12/25 



/ <m 3 (i - t/m 2 D y 2 F D t(b+ -> xe+e~ 

J <m 3 (i - t/m 2 B )- 2 F B T(D° -> Xtv) 



(63) 



where the integrals are over final state phase space, and T is the appropriate integrand 
(which is independent of the decaying meson mass in the approximation that final state 
masses are neglected). The ratio of integrals will depend on the masses of the final state 
mesons: if the integrals are dominated by t ~ t max , then their ratio could be of order 
m 2 D /m 2 B) giving the constraint (j60|) . If instead, the small momentum transfer region 
is dominant, the ratio should be ~ 1. We can skirt these difficulties to some extent by 
comparing B + — > K + £ + £~ to D° — > n + ev because m^jmr) ~ mx/m^. If we assume that 
f+^ K (0) = f+~**(0) [§2J, then we get (|5t]) but with the branching ratio for D° -> ir + ev 
and the upper bound on B + — > K + £ + £~ . For the decays B + — ► n + £ + £~ , we will use (|63|). 
assuming that the ratio of integrals is 1, because this gives the weakest bounds. This 
series of approximations gives the constraints listed in the tables, but we emphasize that 
they are very sloppy. 

CLEO has recently observed the decay B* — > K*^ with a branching ratio of (4.5 ± 
1.5 ± .9) x 10~ 5 ||63|| , which is consistent with the rate expected from the Standard 
Model. Leptoquarks could also contribute via the diagrams of figures 3 and 4, with b 
and s quarks on the external legs, and a charged lepton £ in the loop. We can therefore 
calculate bounds on couplings \ is \ eb by requiring that the leptoquark contribution to 
b — > 57 be "sufficiently small". The leptoquark contribution to b —>■ s'j in a string 
inspired supersymmetric model was calculated in |64j], but is not applicable to the simple 
leptoquark case because of the presence of the supersymmetric partners. 

In section 6.2, we calculate bounds on leptoquarks from the decay fi — > e-y. By 
making the initial fermion a b, the loop fermion a lepton, and dividing by three (no 
colour sum) in equations fl98D, (|101|) , (|104|) and ( |106| ), we can calculate T(b — > 57) at 
a scale // ~ mw ~ fniq from equation (|9"6]). To get a rough idea of the magnitude 
of the bound on leptoquarks from b 57, we can simply add the Standard Model 
and leptoquark amplitudes, run the sum down to m b |StJ, and require that the total 
rate fit into CLEO's allowed window. Since the renormalisation group running mixes 
b — > 57 with b — > sG, we need the leptoquark contribution to b — > sG. Instead, we will 
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assume that it is smaller than the Standard Model amplitude, and neglect it. With these 
approximations, the b — > 57 amplitude at mj can be written 

A(b -> s 7 ) = a[A SM (b -> s 7 ) + .4 LQ (6 -> s 7 ) + 6] (64) 

where a and 6 parametrize the QCD running. The relative sign between the Standard 
Model and leptoquark amplitudes is unknown, as is the exact value of the Standard 
Model amplitude (it depends on m t ), and the relationship between the measured rate for 
B — > K*^j and the calculated one for b — > 37. So to get a rough idea of the magnitude 
of the b — ► sj bound on leptoquarks, we will require 

\ALQ(b^nr,)\&\As M (b^nr,)\ (65) 

where we take m t = 2m\y in the Standard Model amplitude. 
Following reference |65], we get for scalar leptoquarks: 

where £ is a charged lepton, and Qe and Qi q are defined as the quark and leptoquark 
electric charges coming out of the incident b vertex. 
There is a similiar bound for non-gauge vectors: 

^-«<^T£(loitv) 2 (--gauge vectors) (67) 

where 2Qi + \Q lq = 1/3 (2/3) for \ LVo , X RVo and X LVl (Xlv 1/2 , -W 1/2 )- We do not list 
the non-gauge bounds in the tables. The bounds on gauge vectors are to weak to be 
interesting. Equations ( |66"D and fl67l) are really just estimates of the magnitude of the 
bound that could be derived; however, since these couplings are better constrained by 
other B decays (see table 15), this is not a serious problem. 



5 Meson-anti-meson mixing :K°-K°, D°-D° and B°-B° 

Weak interaction box diagrams (see figure 1) neatly account for the small mass term that 
mixes the K° and K° mesons. Similiar box diagrams, with leptoquarks instead of W 
bosons and the internal quarks replaced by leptons would also contribute to the K°-K° 
mass difference, and one can therefore derive bounds on the leptoquark couplings from 
these amplitudes. This has been done in |1(], [L3|, [TJj and more recently in f53|j . 

The Standard Model boxes are GIM suppressed ( of order G ' pa{m q j 'mw) 2 rather 
than Gp<y) because the qqW coupling constant matrix in generation space is unitary. 
This may also be the case for vector leptoquarks, if they are the vector bosons of some 
spontaneously broken gauge symmetry (this is the only renormalisable way of making 
spin-1 particles; however, composite and technicolour models do not neccessarily reduce 
to renormalizable low-energy effective theories.). The gauge boson leptoquark contribu- 
tion to meson-anti-meson mixing would be of the same form as the Standard Model one, 
giving upper bounds on the vector leptoquark couplings of order ([14]) . If the vector lep- 
toquarks are not gauge bosons, GIM-type suppression is unlikely, and the box diagrams 
for vectors would be of order X 4 /mf q . This would give stronger constraints and will be 
briefly discussed at the end of the next section. We assume in the tables however, that 
the vector leptoquarks are gauge bosons, and list the weaker bounds. 

Scalar leptoquark box diagrams will of course be of order A 4 /m^ and we should get 
interesting constraints from these. 
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5.1 Vector leptoquarks 

The Standard Model W boxes contributing to the K°-K° mass difference give [B9 
Am K G F a v 2 m\ 

O tto ^ 7 77 ^ / 7 /7 ^ \ _ _ _ 

+^04^ J ' In -| ] < if | {d-fLs) {d 1(1 Ls)\K > 

m W\ m t~ m c) \ m c/ 

where V is the CKM matrix, and we have used the identity {q^ a ^^ Lq'){q^ ^ ^ a Lq') 
= A{q^^Lq'){q^^Lq'). If one neglects the top quark contributions because Vtd *C 1, and 
approximates the matrix element (vacuum saturation) as 

< K\(d^Ls)(d^Ls)\K >~ l -f K m K (69) 

this agrees quite well with the measured mass difference Amx = 3.5 x 10 -15 GeV. Gauge 
leptoquarks will give the same contribution to Am^-, with gw/V% replaced by Xi q and 
lepton masses instead of quarks. We therefore require 

1 r , x id x is,2™a + x id x is X i'd x i's 2mfmj ^ / mf\ < Gp a sin 2 9 C cos 2 9 C m\ 
32n 2 mf q mf m 2 q (mf — m 2 ) \mj J s/2 47rsin 2 6>vi/ 

(70) 

where I and I' are leptons, and 8 C is the Cabbibbo angle. If we neglect the second term 
on the left hand side, this gives bounds on all vector leptoquarks coupling down-type 
quarks to massive (= charged) leptons I: 

ai\a „ i / m H \ 2 ( l GeV A 



XX < .1 ' g 1^^) (71) 

V100 GeV J V mi ) y 1 

We list the bounds for couplings to muons and taus in the tables. Neglecting the term 
corresponding to the exchange of different flavoured leptons (I and I') in the box should 
be an acceptable approximation, because it would be surprising if it cancelled against 
the first term, and including it makes the bounds more complicated. 

This analysis has been for leptoquarks with left- or right-handed couplings. One 
might hope that for V£ and V^ 2 , who have both, one could escape the GIM-type sup- 
pression by having left- and right-handed couplings at either end of the lepton propagator. 
This would give J2i A^ A^ 0- However, the lepton chiralities must be flipped in the 
propagators to do this, giving a contribution to Am proportional to X 2 L X 2 R mf / mf q , so we 
will neglect these mixed bounds (no stronger; more difficult to calculate). 

There has been no observed mixing between the D° and its CP-conjugate the D°. The 
present upper bound on the mass difference (Amp < 1.3 x 10~ 13 GeV) is considerably 
larger than the Standard Model prediction, but can still be used to constrain leptoquark 
couplings if we approximate the matrix element as in (|69|) (with fo ~ -2 GeV) . This 
gives 



1 lyuycpW < 3Am D 
32tt 2 ' 1 mf Q 2f 2 D m D 



A^ATzx < n£2 u (72) 



implying 

i lu \ Ic 



W\<.s( m 'Q 7— ) (73) 

1 VlOOGeV/ \ m, J K ' 
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A similiar argument for the observed B°-B° mass difference Arris = 3.6 x 10 13 GeV 
gives (using fs ~ -2 GeV) 

m lq \ 2 /lGeV\ 



\ ld X lb \ < .3 



100 GeV ) V 

The limits for muons and taus il = 2, 3 on the internal lines) are listed in the tables. 

We now briefly consider non-gauge leptoquarks. The second term of the massive 
vector propagator: 

9r + ^ (75) 



7 9 2 V XL 11 ' ' 2 

« - m i q rnfg 

is a nuisance in calculations, and the box diagrams we just discussed were done in a gauge 
where it is absent. If the vector leptoquark is not a gauge boson, we do not have this 
option, and the loop momentum integration over the k^k v /mf q terms diverges. However, 
the new physics that produced the vector leptoquark (compositeness?) should appear 
at some energy scale A, which can be used as a cutoff for the divergent integral. The 
non-gauge vector leptoquark contribution to Am is therefore model-dependent. We can 
nonetheless get conservative bounds on the couplings by only considering the contribution 
to the box amplitude from the g^ v terms. This is 

r rl A k k^k u 

Am M = E(A V A i9 )(A^ ) j 



x 



(2tt) 4 {k 2 - m 2 q ) 2 {k 2 - m 2 )ik 2 - m 2 ) 

< MKq^^Pq'Xq^^Pq'^M > (76) 

where M is a K, D or B, q and q' are the meson constituent quarks, P is L or R, 
and i and j are lepton indices. If there are no cancellations in the coupling constant 
sums (no GIM-type suppression), the lepton masses can be neglected and the integral 
approximated as g Ml/ [647r 2 mfJ _1 . Taking i = j, and approximating the matrix element 
as in (H|), we get 

48ir 2 m 2 q 
fhrriM 

for non-gauge vector leptoquarks. This gives the numerical bounds 



. , „ 457T mr 

\\ iq X q | 2 < s2 - lq Am M (77) 



_ 4 / miq \ 
UOOGeV/ 
-a ( m k \ 
UOOGeV/ 
-s ( m k \ 
UOOGeV/ 

where I is an arbitrary lepton flavour index. We do not list the bounds in the tables 
because they only apply if there is no "GIM" suppression. 



X ld x ls 


< 


6 x 10 


x lu \ lc 


< 


2 x 10 


yldylb 


< 


2 x 10 



(KK) 


(78) 


(DD) 


(79) 


IBB) 


(80) 



5.2 Scalar leptoquarks 

As we mentioned, the constraints on scalars are more interesting. The box diagrams 
for a scalar leptoquark (see figure 2) with couplings to leptons of only one chirality (we 
neglect diagrams proportional to (A^A^) 2 ) give a contribution to Am M of 

r d A k k^k u 
Am M = AW* ) / W ( P- S)2( P- m?)( P- mj2) >< 

< M\{q lu Pq'){q^Pq')\M > (81) 
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with the same notation as in (|76|). This is the same integral and expectation value as 
for vectors, but the scalar coupling constant matrix is not required to be unitary, so 
the integral can be approximated as (y' M ^[647r 2 m^ ? ]~ 1 , as in the non-gauge vector case. If 
we take take i = j and assume there are no cancellations between the different lepton 
contributions to Am M , then we get the constraint 

\X iq X iq \ 2 < (82) 

JM m M 



where we have used equation (|69|) to approximate the matrix element. Note that these 
constraints are independant of the lepton mass, so apply equally to neutrinos and all 
generations of charged leptons. We use f D ~ .2 GeV ~ to get the numbers in the 
tables. 



6 Muon Physics 

Muon family number violating processes provide some of the strongest constraints on 
leptoquarks from leptonic physics. Bounds on scalar leptoquarks from muon conversion 
on nuclei, \i — > eee, and ji — > e-y have been considered in H, |^, [II], [IS], K5 \. 

There is a very strong upper bound on the lepton flavour-changing rate for a muon 



to turn into an electron when scattered off a nucleus. The experimental limit [66[ on 
muon conversion on titanium is 

If we neglect the difference in kinematics of the neutrino and the electron, this implies 



< eTi\(fl^Pe)(q llM Pq)\^Ti > 

m lq 



2 <4.6xl0- 12 AGf 



V2 



< uTi\{ji^Liy){u^^Ld)\nTi > 



where P is a chiral projection operator. The matrix elements are not necessarily the 
same, because the isospin structure could be different, and \i-e conversion could be a 
coherent process |)7J]. The ratio of the matrix elements was calculated in [57, for a 



coherent E-/j conversion interaction, and the Nfx Ne matrix element was found to 
be larger than that for TV \i N'v. However, the interaction can only be coherent if it 
is vector, and one can arrange the couplings of some of the leptoquarks such that they 
are only axial vector. We will therefore assume that the matrix elements are comparable 
and require 

-% < 4 x lQ- ll Gl (85) 
m t q 

which applies to all the leptoquark interactions because they all couple u or d quarks to 
charged leptons. As usual, the precise limits are in the tables. 

6.1 fi — ► eee 

Leptoquarks could allow the decay of a muon to eee via box diagrams (see figures 1 
and 2), or through the effective coupling of a muon to an electron and a Z (or a 7*) 
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which subsequently decays to ee. The box amplitudes are similiar to those for meson- 
anti-meson mixing, with quarks and leptons interchanged. The effective Zfie vertex is 
induced by leptoquark triangle diagrams, is divergent, and is non-trivial to calculate. 
By analogy with the Standard Model |70| , we expect the triangle amplitude to be larger 
than the box by a factor of ln(m 2 g /m 2 ); we will nonetheless only discuss the box diagrams 
because they are easier to compute; this will give us conservative limits. 

Scalar and vector leptoquark boxes will induce effective four-lepton vertices 

C LL (fiYLe)(e lfl Le) , C RR ^Re)(e ltl Re) (86) 

In the Standard Model, the muon decays principally (BR~ 100%) to ev e v^ via a W- 
induced four-lepton vertex similiar to (|5(j|): 

^iprfLv^Le) (87) 

If we neglect the electron mass and the difference in the matrix element arising from 
having identical fermions in the final state, we can constrain the effective couplings 
Cll, C rr by requiring 

(Cll,rrY < 8G> 

BR(n — > vve) 

where BR(fi — ► 3e) of course means the experimental upper bound [38| of 10~ 12 . 



.2 ^ 0^2 BR (^ 3e ) 



Had we defined an effective coupling like Cll for the meson-anti-meson mixing, we 
would have had AM/2 = Cll x [hadronic matrix element], so we can estimate Cll and 
Crr from the meson-anti-meson mixing equations. From ( |B"9"D we have for gauge vector 
lepto quarks 



Cll,Crr^-^- 2 (X^X^ (89) 
32tt 2 mf q 

where we have assumed the same quark flavour on both internal fermion lines, multiplied 
by three for colour, and neglected the second term of ([70]) for simplicity. This does 
not give very interesting constraints on the gauge leptoquarks because the amplitude is 
"GIM" suppressed, giving a decay rate proportional to X 8 m q /mf q . However, it does give 
bounds on couplings of electrons and muons to quarks of all flavours (except t). We 
have, from 



where q is a quark flavour other than top. 

For non-gauge vectors we have, in the same approximation used for meson-anti-meson 
mixing: 

3(A^) 3 A 2 ^ 



Cll, Crr — — — (91) 
32vr 2 m 2 g 



This gives a more interesting bound than 

y/X^iX 1 ") 3 < 6 x 10~ 3 ( 1Q( ^Q eV ) ( non - £ au g e vectors) (92) 
where again q ^ t. We do not list the non-gauge results in the tables. 
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Scalar leptoquarks give 

_3(A^ 

12g7r2mj , (93) 

for light quarks, which implies (g 7^ i) 

y/\*>W < 10- 2 (y^^) (scalars). (94) 

These approximations for C^x and (Trr only apply if the internal line fermions are 
light. For meson- anti- meson mixing induced by leptoquarks, this is a good approxima- 
tion because the internal line fermions are leptons. However, one could in principle get 
constraints on leptoquark couplings to the top quark from fi — > eee box diagrams involv- 
ing top quarks and leptoquarks. The exact expression for the intergral in ( |76| ) and ( PD 

1 ' ' ' ' + n ^ \2 ^ ax v~ 7i~~^~T2 ^ n x t 1 (95) 



647r 2 m 2 g x 4 — x g 1 — x g 1 — x t (1 — x q ) 2 9 (1 — x t ) 2 

where x t = m 2 /mf q , x q = m 2 /m 2 q , and q is any up-type quark. Any bounds will clearly 
depend sensitively on the top and leptoquark masses, so we neglect them here. 

6.2 \i — > e7 

The lepton flavour violating decay /i — > could be mediated by leptoquark triangle 
diagrams (see figures 3 and 4). This decay will again be "GIM" suppressed for gauge 
leptoquarks, but in this case will be of order \ 4 m q /mf q rather than 0(X 8 m q /mf q ) as in 
the n — > eee case. There are also fewer powers of A in the scalar and non-gauge vector 
rates. 

The decay rate of an initial fermion into a photon and a final fermion ff is [[fl] 

r(/< - fn) = ? (1 - — ) 2 (1 - (— ) 2 ) [^ y (o) 2 + ^ A (o) 2 ] (96) 

where F v and F A are the vector and axial vector coefficients of the magnetic moment 
term in the matrix element: 

{n^ v n 
... + —^{F v (q 2 ) + F A ( q 2 ) l5 } + 
rrif + vrii 

and a^ u = ^[7^,71- 

The \ 2 L and \ 2 R terms in the scalar leptoquark contribution to (g — 2) e have been 
calculated in detail in 0, [75|, Djouadi, Kohler, Spira and Tuta [f76| , and Morris |74 




also calculate the contribution proportional to Xl^r- These should be closely related 
to Fy and Fa] if one set rrii = rrif = m e in the final result for Fy, one ought to get 
e(g — 2) e = Fy. Unfortunately, in this limit our approximate result is smaller than the 
exact one in [[T6| by a factor of 1/3, and has a different coefficient in front of the leptoquark 
charge. In spite of this disturbing discrepancy, we will use our formula because it gives 
weaker bounds on the couplings, and because it is three times the result in |75| (three 
for colour), whose results can be checked against Standard Model calculations. (Davies 
and He jff5f seem to multiply Leveille's result by 9 in their published paper, which would 
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agree with J76]. This is apparently a typographical error, and should have been a three 
|77| .) Note that we have also neglected QCD corrections to the diagrams. 

If we assume m q <C rrii q (this assumption was not made in J7B], 75 1, so their result is 
applicable to all quarks, in particular the top), we can approximate 



A! 



rrii + mf rrii + rrif (87r) ; 



mi±m f Qi q 
-ei—r 1 - Q q ) 



(scalars) 



(98) 



for couplings of the same chirality at either end of the scalar leptoquark propagator. (The 
same bound applies to \ 2 R .) Q q (Qi q ) is the quark (leptoquark) electric charge coming 
out of the incident muon vertex, "e" is the electromagnetic coupling constant, and X lq is 
an orthogonal matrix in our approximation that all coupling constants are real. The + 
and — on the RHS of ( |9"8"D are respectively for F v and F A (hence the subscipts on F v 
and F A ). 

The X^X R contribution (from scalar leptoquarks) can be larger than the X 2 L , X 2 R 
terms, because it is proportional to (mim q /m 2 q ) \n[m 2 q /m 2 q }. We will nonetheless neglect 
it because it only applies to those leptoquarks which have couplings of both chiralities. 



These terms are present in the calculation of 173, |76 



Requiring that radiative muon decay, as calculated from (|96|) and (|98|) take place 
more slowly than the present experimental upper bound r(// — > ej) < 1.5 x 1CT 29 GeV, 



gives 



X 2 T q X^, X 2 R q X R < 8 x 10- 5 ( r^-^Y 



(scalars) 



(99) 



VlOOGeV/ 

where q is a quark other than the top, and we have assumed here (but not in the tables) 
that ~l — Q q ~ 1. For the SU{2) doublet and triplet leptoquarks, this formula is an 
underestimate, because it assumes that only one leptoquark is being exchanged (rather 
than two or three of different charges). For simplicity, we estimate the rate as being due 
to the member of the mutiplet that gives the largest (^ — Q q ). For £1/2, (3 — Qi) = ® 
(remember that Q q and Q\ q are defined as the charges coming out of the first vertex), 
so we have no bound on the tables. There will certainly be a contribution to Fy from 



S1/2 leptoquarks, but it will be suppressed with respect to fl98|) by powers of m q /mi q or 
m n/ m iq, so we neglect this. (It is also possible that ou approximate formula is wrong, 
and F v ~ aQi q + Q q with a 7^ —1/2, in which case the bound would be of order (|99|).) 
The two triangle diagrams for gauge vector leptoquarks contribute approximately 



X f L q X l le (m, ±m f ) 



mi + rrif mi + rrif 
2 



m 2 7 m n 
-f(- + 31n-f 
mf q 2 mf q 



+ 



327T 2 



13 m\ 
4 m\ 



x 



m 



lq 



(gauge vectors) 



(100) 



(101) 



where we have neglected a number of (presumably smaller) terms, and the same formula 
would apply to couplings X%. The term of order [XX T ]if /mf q vanishes fori^f due to 
the orthogonality of X lq . 

we can check this expression by comparing it to Leveille's 



If we set rrii 



rrif 



m 



[ 75| formula for the massive gauge boson contribution to g — 2. Requiring 



r(// -> ej) 



a 



L 



in 



(4tt) 



m 



lq 



q -K 



Qq 



2 mfq 



13 



Qlq 



< 1.5 x 10"' y GeV 
(102) 
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gives 



\2<j<vlq . -2 / Wllq \ A ( 1 GeV 

A L A L i A R A R < 



VlOO GeV 



m 9 



(gauge vectors) 



(103) 

for gauge vector couplings. This applies to all quark flavours except top. Due to the 
number of approximations in this calculation, this constraint can only be considered a 
rough estimate. 

It would not be difficult to modify the Standard Model amplitude for b — > S7 to 
descibe a leptoquark and a quark of arbitrary charge in the loop. Such a "modified 
Standard Model" expression could be used to constrain gauge vector leptoquark couplings 
to the top quark. However, as in the case of meson anti-meson mixing, the formula is 
fairly complicated, and does not provide clear constraints on the coupling constant when 
both the top and the leptoquark masses are unknown. 

Some of the vector leptoquarks (V^ 1 , Vy 2 ) have two different interactions with Stan- 
dard Model fermions. The two couplings have opposite chiral projection operators. We 
argued that for the box diagrams, the bounds on the product XlXr would be of the same 
magnitude as the "GIM" suppressed bounds on A^A^ or XrXr because the fermion mass 
must appear in the amplitude to flip the chirality between the gauge boson vertices. 
Since the effective coupling for the four- fermion vertex has dimension mass -2 , it must 
be of order X^X R mj/mf q , where mj is the mass of the internal line fermions. This is 
roughly the same as the "GIM" suppressed amplitude constraining X\ or X 2 R . However, 
this is not true for the triangle graphs, where the bounds on A^A^ are much better than 
those on A| or X R : the coefficients F v and Fa are dimensionless, so putting a quark mass 
in the numerator to flip the chirality simply changes the leading contribution to Fy and 
Fa from X\rr? li jm l q to XLXRm^m q /m 2 q . Since there is no reason to expect Y^ q X^X e R = 



for vector leptoquarks, we have J78 



rrii + mf rrii + mf IQtc 2 mf c 



lj( n ' q (X l ^±X^X^(Q lq -Q q ) (vectors) (104) 



lq 



where e is the electromagnetic coupling constant, and i and / respectively label the 
initial and final leptons. Taking (Qi q — Q q ) ~ 1 (we do not assume this in the tables), 
this gives 



AM, AgA? < 10- 6 ( 2 ( 1 -^) (vectors) ( ] 00 ) 



100 GeV/ V m i 



from the upper bound on the decay /i — > ej. 

For non-gauge vectors, the leptoquark propagator is again of the form (75), and the 



second term will give cutoff dependant contributions to the decay rate. If the leptoquark 
mass is well below the scale at which new physics appears, these terms could be large. 
We will nonetheless neglect them, as we did for the boxes. This should make our bounds 
conservative. So using only the g^ v part of the vector propagator, and assuming that 
there are no cancellations between triangles containing different quark flavours, we get 
the leading order terms of ( 101p : 



+ 1 \ F -\ We (mi±m f \ , n „ 5^ . , 

- - 1 J -)(2Q q + -Q lq ) (non - gauge vectors). 



rrii + mf rrii + mf 32ir 2 \ m 



(106) 
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This gives 

\ 2q \ lq < 2 x KT 5 ( — — \ (non - gauge vectors) (107) 
V100 GeV/ 

where we have again taken (2Q q + \Qi q ) ~ 2. 

7 r decays 

7.1 Semi-leptonic decays 

The tau is the only lepton heavy enough to decay to a meson and another lepton. The 
Standard Model lepton family number conservation only allows decays to v T + X, and 
t — > ixv is observed at the expected rate. The leptoquark couplings Xls , ^LSi, ^lv an d 
\lvii which couple a tau, a neutrino, and an up and down type quark, must therefore 
satisfy 

e 2 AG F , . 

< < 7T • (108 > 

Leptoquarks could also mediate r — > eM, fiM, where M is a neutral meson lighter than 
the tau. If we assume that 

< 0\u i 'f'y B u i \M o >~< Ol^j^d^M >~< Olu^V^I-W" > ( 109 ) 

(for instance, < 0|w7 M 75w|7r° >~< 0|d7 At 75u|7r + >), and neglect the mass of the de- 
cay lepton, then we do not need to calculate the tau decay rate to a lepton (£) and a 
meson(M). For a four- vertex of the form 



e 3k e ij 



(f 7 W)(?V9 fc ) (HO) 



we can simply require 



v m lq J BR ( T M ^ ) 



where V is the appropriate CKM matrix element for the decay t~ — > M~u, k and j label 
the flavours of the constituent quarks in M, and BR(r~ — ► M°t~) is the upper bound on 
this process [73 ]. Using t = e, fi and M° = 7r°, .fT , we get the constraints listed in the 



tables. Note that unlike the pseudoscalar meson decays, there is no chiral suppression so 
we would not get better bounds on the couplings of the effective pseudoscalar vertices. We 
therefore do not compute constraints on the product of left- and right-handed couplings. 

7.2 Leptonic decays 

Leptoquarks will contribute to the decays r — ► III and r — ► £j {£ — fi, e) in the same 
way as they did for muons. The bounds however, will be weaker, because the rare decays 
of the tau are not as tightly constrained as those of the muon. Using equations (p3|). 
( j9lD and BR(r -> 3/) ^ 1.5 x 10~ 5 {7pJ, we find that for scalar and non-gauge vector 
leptoquarks 

A?(A^) 3 , A^(A^) 3 < .2 (y^I^) 2 (non - gauge vectors) (112) 
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AM) 3 , A^) 3 < .3 (scalar) (113) 

where q is any quark flavour other than top, and I can be \x and/or e. The bounds on 
gauge vector leptoquarks are to weak to be meaningful. 

The experimental upper bounds on the decays r — > £7 are T(r — > ^7) < 9.2 x 10~ 18 
GeV |72[ and T(r — > e7) < 2.6 x 1CT 16 GeV ||73[. The radiative muon decay bounds 



can be easily modified to constrain leptoquarks coupling to the r. The constraints on 
gauge vectors with couplings of one chirality are again to weak to be worth listing, but 



equations (|96|) and ( |104j ) give bounds on vector couplings of both chiralities (Xl^r) 



which are in the tables. One can also get bounds on non-gauge vectors from (|106|) and 
the upper bounds on r — > pq and r — > ey. This gives 



a?a?, a?a* < Q^n^ (io^v) ( non - ^ e vectors ) ( 114 ) 

A?A?, AjAj < ( T ^ v ) 2 (non - gauge vectors) (115) 

Equations (|96])and (p8|) can be used to constrain scalar leptoquarks contributing to ra- 
diative tau decays. We list the bounds in the tables. 



8 Other Constraints 

8.1 Quark- lepton universality 

As was pointed out in 0, leptoquarks would contribute at tree level to neutron j3- 
decay, but only via box diagrams to // — »■ euv. We know from (|89"D , (|91|) and (|9~3"D that the 
effective leptoquark induced coupling constant for the 4-lepton vertex is much smaller 
than that for the 2-quark 2-lepton vertex: 

C LLtRR « ^ (116) 

m lq 

We can therefore neglect the leptoquark contribution to muon decay, and obtain con- 
straints on the udev couplings from requiring that the Fermi constant in f3 decay not 
differ signifcantly from the muon decay measurement. If the effective Fermi constant for 
P decay =Gp was 

^ = ^ + 4 cur) 

V2 y/2 m% 

then the rate for n —>■ peu would be increased by a multiplicative factor 



p 2 p 4 
+ V2G F m% + 8Glmf q 1 ) 

The CKM matrix element V ud is determined to be 0.9744 ±.0010 @ by comparing Gp 
to Gp, and this value is consistent with unitarity. We can therfore require that the 
leptoquarks contribute less than the error in V u d- 



V2e 2 L 5V ud 
Am? K 



< -y^G f (119) 



Hq V ud 
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where we take 5V u d = 0.9754— [the minimum value given in the Particle Data Book] 
|1. This gives 

for udev axial vector couplings. We neglect the contributions from effective scalar ± 
pseudoscalar quark matrix elements (~ XlXr), because they are strongly constrained by 
the ratio R v . (They have also been carefully considered in 



8.2 g - 2 

Since the magnetic moment of the electron and the muon are two of the most accurately 
measured quantities in physics, one always hopes to get good constraints on new particles 



from them. Constraints on leptoquarks from g — 2 have been calculated in |LJ| and |76 
It is easy to see from equation (|97|) that 

(9 - 2 )e = m e F v 
2 e rrii + nif 

with rrii = rrif = m e . Using Fy/(mi +rrif) from equations (|98]), ( |101| ), ( |104j ), and ( |106| ). 
we can constrain some of the couplings by requiring that the leptoquark contribution to 
g — 2 for the muon and the electron be less than the difference between the theoretical 
predictions and the experimental measurements. Although [g — 2] e is more accurately 
measured than [g - 2]„ (A[g - 2] e /2 = 6 x KT 10 , A[g - 2\J2 = 2 x 1(T 8 ), the lepto- 
quark contributions are proportional to the lepton mass, or mass squared, so the relative 
smallness of m e suppresses most leptoquark contributions to [g — 2] e . As can be seen 
from the tables, these bounds are considerably weaker than those from radiative decays, 
but constrain different combinations of generation indices. We have neglected couplings 
to the top, because we assumed m q <C rrii q in calculating Fy. As previously mentioned, 
this approximation was not made by Davies and He [15|, or Djouadi, Kohler, Spira and 
Tutas fT6|| , who give a full analytic expression for the scalar leptoquark contribution to 
g — 2, and a reference [79|] for the non- gauge vector contribution where the divergent 



k^k u /mf q term in the propagator is included. 

From equation ( |101| ), we can get weak constraints on the gauge vector leptoquarks 
that contribute to g — 2. Since the initial and final lepton are of the same type in this 
diagram, there is no "GIM" suppression, and the leading order terms in ( |101| ) give the 
bound 

^«<^(lo^) 2 <-) 

where q is any quark other than top. The bounds on A^Ar for gauge vectors are stronger 
than this (for heavy quarks in the loop). We get, from (g — 2) e and (g — 2)^: 

A " A? K Q^Q, ( W&v) (123) 
AM< Q^(l^v) 2 ■ (124) 

The bounds on scalar leptoquark couplings of the same chirality (A|, A^) from (g — 2)^ 
are similiar to ( |122| ). As in the vector case, we could probably get stronger bounds on 
Xl^r, at least for heavy quarks on the internal lines. We neglect this, because there are 
better constraints from elsewhere (see table 15), and bounds on the combination A^A# 
are less interesting in the first place. 
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8.3 Neutral current couplings 

The many different experimental determinations of sin 2 9\y seen to agree quite well. 
Since leptoquarks mediate two-lepton, two-quark interactions like the weak gauge bosons 
of the Standard Model, they could contribute at tree level in some, but not all, of these 
experiments. There will therefore be bounds on leptoquark couplings from requiring that 
they not disrupt the agreement between various determinations of sin 2 9w- These will be 
particularity interesting, because they apply to leptoquarks coupling to first generation 
fermions- precisely the leptoquarks that could be seen at HERA. 

Leptoquarks should not have a significant effect on the LEP measurements, because 
LEP runs on the Z mass. However, they should contribute on an equal footing with the 
electroweak gauge bosons in atomic parity violation experiments, and in deep inelastic 
scattering of neutrinos off nucleons. It turns out that the best bounds on the largest 
number of leptoquark coupling constants come from atomic parity violation experiments 
on cesiumj38|, and these constraints have been considered by Langacker | 8~3|| . The atomic 
weak charge Qw = —2[Ci u {2Z + N) + Cid(Z + 2N)] of the cesium atom is measured to be 
Qw = —71.04 ± 1.58 ± .88 (the second error is theoretical). C\ u and C\d are coefficients 
in the effective four-fermion parity-violating Lagrangian 

G 

-^{h fi j 5 l){Ci u uj fi u + C ld dj^d) (125) 

where at tree level (see |3j| for the one-loop Standard Model corrections) C\ q = —I 3 + 
2Q em sin 2 6w + AC\ q) and AC\ q is the leptoquark contribution: 

e 2 \/2 

m lq 4Gr F 

The Standard Model expectation for Q w is —73.21 ± .08 ± .03, so requiring that the 
leptoquark contribution to Qw be less than the difference between the experimental and 
Standard Model values, gives 

|ACi tt | = |AC W | < .003 (127) 

for leptoquarks that induce AC\ U = ACid, and 

|Ad«|, |AC ld | <.006 (128) 

for leptoquarks where ACi u = or AC^ = 0. These bounds are slightly weaker than 
those in |83||. They also neglect the sign of the corrections, which is conservative, but 



not altogether reasonable because the theoretical prediction is one and half a to one 
side of the experimental result. From (|127|) , ( |128|) and tables 3 and 4, one gets that for 
Alv i; ^rv 1/2 , ^rs 1/2 and Als x 

e2< - 001 (l^v) 2 (129) 

and for all the other couplings that only contribute to AC± U or to ACid- 

■ (130) 

lOOGeV J V ; 

We have assumed for these constraints that X RSo ^ X LSo , Xrs 1/2 ^ls 1/2 , >^rv ^ >^lv , 
and XlVim 7^ ^RVi/2i because the parity violating coupling is A^ — A#. This is a reasonable 



assumption, as noted in S3 , because there are strong bounds on the product of these 



couplings XlXr from the ratio (it — > eu)/(ir — > (see section 4.1) 
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8.4 Neutrino Oscillation Experiments 



Four of the leptoquark couplings induce four-fermion interactions that couple a neutrino, 
a charged lepton and an up and down type quark (see tables 3 and 4). The couplings 
\ls i ^lv an d A^Vi can therefore be constrained by neutrino oscillation experiments 
that look for charged leptons of a different flavour from the original neutrino beam. 

If the neutrino masses are small, and one neglects the mixing of vi a and v\ h with the 
third neutrino, then the probability that a neutrino v\ a with energy E will oscillate into 
a vi h over a distance x is 

|<^|^>| 2 = sin 2 2^sin 2 (^) (131) 

where 9 is the angle mixing z// a and ui b , and Am 2 = \rn 1 Va — m. 2 J. If sin 2 (Am. 2 a;/4.E) 
oscillates rapidly enough, (large Am 2 ), it can be averaged to 1/2. The "large Am 2 " 
bounds on the mixing angle sin 2 29 are therefore bounds on the probability for a v\ a to 
turn into a Ui b , or equivalently on the ratio of l a to l b produced by the neutrino beam. 
This can be used to constrain the leptoquark induced rate for v\ a + iV — > l b + X {pi a +u — > 
l b + [d, s, b}) and v Xa + N — »• l h + X {v ia + u -> l b + [u, c]): 

e 2 

— \ < 2G F \sm29\ (132) 
m i q 



where sin 29 is from the "large mass" upper bound from accelerator experiments |38j . 

The scalar leptoquarks that interact with neutrinos have fermion number 2, so couple 
an incoming v to an incoming d, or an incoming v to an outgoing d. The vectors, on the 
other hand, have fermion number 0, and couple an incoming v to an outgoing u or an 
incoming v to an incoming u. This is important because we assume that the outgoing 
quark can be any flavour other than top (we neglect the phase space suppression due 
to the mass of the b), so to get the flavour indices right on the couplings constants, we 
need to know if the neutrino couples to the leptoquark and the outgoing quark, or the 
leptoquark and the incoming quark. 

The bounds are listed in the tables. 



8.5 Asymmetries at colliders 

The constraints on leptoquarks from the forward-backward asymmetry (Apg) at ex- 
isting ee colliders (excluding LEP1) were considered by Hewett and Rizzo JJ7). (Lepto- 
quarks are not constrained by Afb at LEP1 because it runs on the Z mass.) Possible 



constraints from future colliders were considered in |85 1 . Hewett and Rizzo required that 
Afb and the total cross-section for ee — > qq not differ from the Standard Model values 
by more than 5 or 10% at ^fs = 40 GeV. They claim that a greater contribution would 
disrupt the agreement between Standard Model predictions and the measured rates for 
b and c quark production at PEP and PETRA. 

Their bounds apply to the couplings (eu)(eu) and (ed)(ed) ( (11) (11) in the tables), 
for scalar leptoquarks. A similar constraint could be calculated for vector leptoquarks, 
but since the measurement of sin 2 9w in parity violation experiments gives more stringent 
limits, this is unneccessary. We list Hewett and Rizzo's bounds in the tables as 



2 



A2 < 02 (isiv)' • < 133 > 



31 



Note that this analytic bound is just an approximation to the limit graphed in their 
paper. 

In principle, some leptoquarks could contribute strongly to Alr at tree level, because 
they only couple to electrons of one chirality. However, measuring Alr requires polarized 
electron beams, so we can not get bounds from this at present (SLC has polarized beams, 
but, like LEP1, is running on the Z mass). 



8.6 Muon beams 

As was noted by Heusch one should in principle be able to get interesting bounds 
on leptoquarks from muon beam experiments. Leptoquarks contributing to fiN — > \iX 
can only be constrained by measurements of asymmetries in this process, because the 
rate is dominated by photon exchange. However, if one could effectively search for 



/iiV — > £X |81 |, where £ is an electron or a tau, one could constrain couplings with the 
family structure \^ q \ lq , where q and q' can be any type of quark (other than top) if one 
allows the muon to scatter off sea quarks in the nucleon. 

We can get weak constraints on leptoquarks mediating the reaction fiq — > \iq' from 



the measurement [36] of the asymmetry 



B = ~ d f[j\ (134) 

where da(n] i ) is the cross-section for fiN —>■ fiX for an incident right-handed muon. If 
one assumes that only valence quarks contribute to this process, one can show that the 
interference between Z and 7 exchange gives |86|. |H7| 



B = -PCy)(c d A - 2c\f-g(y) (135) 

where Q 2 is the energy transfer, p is the polarization (= .81 in this experiment), g(y) 
is a function of the ratio of outgoing to incident muon energies (y = 1 — Ef/Ei), and 
Cy , are the usual vector and axial vector couplings, normalized such that = —1/2. 
This formula can easily be modified to descibe the interference between leptoquark and 
photon exchange by substituting e 2 L /m 2 q <-» A^/2Gp, and setting the various c's to ±1 
or as required. We then require the leptoquark-photon intererence contribution to the 
asymmetry (|135|) to be less than the error in B. This gives the bounds listed in the 
tables. 



9 Discussion 

We have attempted to make this catalogue of constraints as comprehensive as possible. 
However, it is inevitable that we will have omitted some (we have, for instance, neglected 
bounds from CP violation). One could also do the calculations more carefully than we 
have done, and probably improve many of the constraints, as we have tried to err on the 
side of caution. 

In the tables, we list the constraints on gauge boson vector leptoquarks. In some cases, 
the bounds on non-gauge vectors are much stronger, and we have tried to estimate these 
in the text. 
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The data in tables 5 to 14 is difficult to interpret, so we have collected the best 
bound on each product of couplings, and the experiment it came from, in table 15. The 
couplings Xls , ■■■^lv 1 from equations (1) and (2) are really 3x3 matrices in generation 
space. We therefore list the matrices horizontally and the indices of the two matrix 
elements whose product is constrained vertically. So the top right-hand entry of table 15 
is the best bound on A" s Aj^ . As before, the lepton index comes first, and the numbers 
in the table are multiplied by (m; 9 /100 GeV) 2 . We have assumed, for this table, that 
all neutrinos are light (m v <C m e ). ALthough we calculated these limits for leptoquarks, 
many of the bounds would apply to other massive "beyond the Standard Model" bosons. 
In particular, the rare meson decay bounds apply to any particle that induces L V ± A 
lepton current' x ' V ± A quark current' four-fermion operators. 

Ideally we would like to have individual bounds on the squares of each element of 
every matrix. Instead, we find that the squares of each matrix element are in general 
unconstrained. (By this we mean that for a leptoquark of mass mi q ~ 100 GeV, A 2 
could be of order 1.) This is hardly surprising, since leptoquark-mediated interactions 
whose matrix elements involve the combination X lq X lq do not change lepton generation or 
quark flavour; with two exceptions, they are therefore competing with electromagnetic 
amplitudes for the same process. For instance, leptoquarks could mediate the decay of 
qq mesons, but since the photon can also, there are no interesting bounds on leptoquarks 
from these decays. The two exceptions are parity violation experiments, which is where 
the bounds on (A 11 ) 2 come from, and if one of the leptons is a neutrino. Al5 d , Xlsi, Alv d 
and Xlv! are the only couplings to involve a neutrino, a charged lepton and an up and 
down type quark, and are sometimes better constrained in consequence. 

Most of out constraints come from processes that are suppressed or forbidden in the 
Standard Model by lepton family number conservation, or the absence of flavour-changing 
neutral currents among the quarks. Such limits tend to involve the product, rather than 
the square, of couplings: A^A mri < something, with (ij) ^ {ran). An upper bound on 
the product of two couplings is not terribly useful, because either coupling could be very 
large or very small. The other problem associated with bounds from flavour changing 
processes is that usually i ^ j and/or m ^ n in bounds of the form A* J A mn < something. 

If the leptoquark couplings have a generation structure similiar to the Standard 
Model, we might expect couplings \^ with i ^ j to be one or two orders of magnitude 
smaller than the i = j value, by analogy with the CKM matrix. It is unfortunate that the 
couplings that we might expect to be largest are the ones we have the fewest constraints 
on. 

We have few bounds on leptoquark couplings to the top quark, because these can 
only be obtained from loop diagrams (meson-anti-meson mixing and /i and r decay to 
three leptons or a lepton and a photon), and the quark mass was neglected in the triangle 
loop calculations. This could be corrected, but would still not lead to simple numerical 
bounds, because they would depend on the unknown leptoquark and top masses. The 
interactions which are unconstrained because they involve a top quark are indicated by 
stars in table 15. 

It is worth noting that some of the interactions couple a charged lepton to an up- 
type quark, some to a down, and some to both. The limits on the coupling constants of 
interactions involving down-type quarks tend to be better than those for ups, because 
Ks are better constrained than Ds, and Bs exist. 

HERA expects to be able to see leptoquarks with mi q ^ 300 GeV and A 2 ^ 10~ 4 
as peaks in the e ± p — > e ± X cross-section as a function of x 0, so if the upper bound 
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on A* 1 A™ 71 from a rare decay is greater than 10 _5 (m/q/100 GeV) 2 , then the leptoquark 
'a' can be produced and detected via these couplings at HERA. As is clear from table 
15, this means that it is consistent with the bounds computed here for HERA to see 
any of the leptoquarks. Four of the constraints listed here are strong enough to each 
rule out a particular production and decay channel, but since the leptoquark could 
easily decay to a different quark-lepton final state, this is not very significant. For 
instance, the bounds from muon conversion on titanium (e 2 < 7 x 10 _7 (m^/100 GeV) 2 ) 
suggest that no leptoquark produced in the collision of an electron (or e) and a first 
generation quark will be detected decaying to a muon and a first generation quark (a 
second or third generation quark would, however, be quite possible). Similiarly, the 
absence of kaon decays to p,[i, nfie or ep, implies that HERA will not see the leptoquarks 
VoS ^1/2 > ^i > S\/2, and Si decaying to a strange quark and an electron or muon (but 
again, some other flavour of quark would not be ruled out). If the leptoquark Yukawa 
couplings were generation independant, these constraints would apply to all initial and 
final quark-lepton combinations, and HERA would not see leptoquarks. However, if one 
makes the far more reasonable assumption that the couplings are generation dependent 
(very true of the Standard Model fermion-fermion-boson couplings), then our bounds 
from rare processes do not seriously infringe on HERA's chances of seeing leptoquarks. 
The experiments listed above do constrain e 2 /mf q < 10~ 9 GeV~ 2 (e 2 < 10~ 5 in the 
tables), but they do not apply to the couplings for e ± p — > e ± X (where X consists of u 
and d quarks). So leptoquarks can easily have masses and couplings that are accessible 
to HERA and consistent with the upper bounds on rare processes. 

In figure 5 we have plotted, following and using quark distributions from J88| , 
the differential cross-section for e~p — > e~X in the presence of the scalar leptoquark 
S Q with a mass of 200 GeV, and couplings A^ = .008, A^ 5o = for i,j ^ 1, and 
Ars d = 0. This choice of mass and coupling is consistent with the upper bound (A^) 2 < 
2 x 10 _3 (m/ (? /100 GeV) 2 from quark-lepton universality. One could make similiar plots 
for other leptoquarks, and second or third generation final state quarks and leptons. It 
is clear that if ZEUS and HI do not see leptoquarks, they will impose strong constraints 
on all couplings of the form A 11 A mn . 



10 Conclusion 

We have calculated generation dependent upper bounds on the coupling constants of all 
renormalizable, B and L conserving interactions consistent with the SU(3) x SU(2) x U(l) 
symmetry of the Standard Model, involving a lepton, a quark, and a scalar or vector 
leptoquark. The constraints come from rare meson decays, meson-anti-meson mixing, 
lepton decays, and a miscellany of electroweak tests. We list the bounds on the square of 
each coupling constant from each experiment in tables 5 through 14; since the constraints 
depend on the leptoquark mass, and only apply for certain generation indices, we quote 
numerical bounds based on mi q = 100 GeV, and put the lepton and quark generation 
indices in parentheses ( so that the top right hand corner of table 5: (11) (nl) < 10~ 3 
means that Af^A-Qj < 10 _3 (m ig /100GeV) 2 ). In table 15, we list the best bound on each 
coupling constant for every combination of generation indices. 
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